Magnetic ordering and spin waves in Na .82CoO2 
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Na^Co02, the parent compound of the recently synthesized superconductor Na^CoC^yH^O, 
exhibits bulk antiferromagnetic order below ~ 20 K for 0.75 < x < 0.9. We have performed 
neutron scattering experiments in which we observed Bragg reflections corresponding to A-type 
antiferromagnetic order in a Nao.82CoC>2 single crystal and characterized the corresponding spin- 
wave dispersions. The spin waves exhibit a strongly energy-dependent linewidth. The in-plane and 
out-of-plane exchange constants resulting from a fit to a nearest-neighbor Heisenberg model are 
similar in magnitude, which is unexpected in view of the layered crystal structure of Na^CoC^. 
Possible implications of these observations are discussed. 

PACS numbers: 75.30.-m, 76.75,+i, 72.80.Ga, 71.30.+h 



The cobaltate Na^OoC^ : 2/H2O has recently enjoyed 
intense attention. The composition with x ~ 0.30, 
y ~ 1.4 has been shown to be superconducting over a 
narrow range of x, with maximum transition temperature 
T c ~ 5 K 1, 2]. This compound is particularly interest- 
ing because its structure is similar to that of the high-T c 
copper oxide superconductors. In both materials, super- 
conducting sheets containing oxygen and a spin-1/2 tran- 
sition metal are separated by layers of lower conductivity 
in an anisotropic crystal structure. However, a number 
of characteristics suggest that the superconductivity in 
this compound may be unusual in different ways from 
that found in the cuprates. For example, some exper- 
iments indicate that the symmetry of the Cooper pair 
wavefunction may be p-wave 0, Q . 

The unhydrated parent compound Na K Co02 is inter- 
esting in its own right owing to its exceptionally high 
thermopower over the range 0.5 < x < 0.9, which, un- 
usually, accompanies low resistivity and low thermal con- 
ductivity []| |g. Magnetic susceptibility data for com- 
pounds with 0.5 < x < 0.7 show Curie-Weiss behavior 
with a negative Weiss temperature 0, @, El • No static 
magnetic ordering has been observed for x < 0.7 in /iSR 
experiments ^(J . The magnetism in this system depends 
sensitively on the doping level: /iSR experiments per- 
formed on powder samples with x = 0.75 suggested the 
presence of a magnetically ordered phase below T = 22 K 
nil. Recent magnetic susceptibility measurements have 
shown evidence of AF long-range order below 20 K for 
0.75 < x < 0.90 [EII3 Anisotropic DC magnetic suscep- 
tibility and /iSR measurements on the x = 0.82 composi- 
tion both showed that the Co spins are oriented along the 
c-axis [12|. Searches using unpolarized neutrons for ev- 
idence of corresponding static magnetic order have thus 
far been unsuccessful [l3j . 

In a recent time-of- flight experiment on a Na.j;Co02 



crystal with x — 0.75, Boothroyd et al. observed fer- 
romagnetic (FM) fluctuations within the afe-planes [l3| . 
However, the c-axis momentum transfer could not be var- 
ied independently of the energy transfer in this experi- 
ment, and correspondingly only a projection of the scat- 
tering cross-section onto the a6-plane was probed. The 
magnetic ordering pattern below Tn, as well as the ex- 
change parameters and their relationship to the macro- 
scopic susceptibility above Tn, have thus far remained 
undetermined. These issues are addressed in the neu- 
tron scattering studies reported here. For single crys- 
tals with a;=0.82, we find static ordering and low-energy 
fluctuations characteristic of A-type antiferromagnetism: 
namely, antiferromagnetically coupled ferromagnetic lay- 
ers. Surprisingly, the AF interlayer exchange constant 
resulting from a fit to a simple Heisenberg model is al- 
most as large as the intralayer FM exchange, despite the 
difference in Co-Co distances. A possible relationship of 
this observation to the Co 3+ -Co 4+ charge-ordered states 
recently suggested on the basis of NMR 0, 0] and op- 
tical conductivity jl6l | experiments is discussed. 

Single crystals of 7-phase Naa;Co02 were grow n by the 
floating-zone technique in an image furnace |l2l Il7| . The 
sodium and cobalt contents were determined through 
analysis of pieces cut from the same ingots and ad- 
jacent to the crystals used for the neutron scattering 
experiments, using inductively coupled plasma atomic 
emission spectroscopy (ICP-AES) and atomic absorption 
spectroscopy (AAS). The Na:Co ratio was found to be 
x = 0.82 ± 0.04 for each of the two crystals examined in 
this experiment. Magnetic susceptibility measurements 
were performed on portions of the ingots adjacent to the 
crystals used in the experiment. All of the pieces tested 
exhibited the expected AF transition at ~20 K [l^ . 

We discuss first the determination of the magnetic or- 
dering pattern by elastic neutron scattering. We used 
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FIG. 1: Normalized SF intensity (equal to the SF intensity 
divided by the NSF intensity, viz. the inverse of the flipping 
ratio as it is usually defined; see text) at Q = (101) and (100), 
as a function of temperature. The data were taken with the 
neutron polarization P || Q, with ki=2.QQ2 A -1 . The lines are 
a guide for the eye. Inset: the A-type AF structure, shown 
with Co spins || c. 



the polarized-beam spectrometer 4F1 at the Laboratoire 
Leon Brillouin, Saclay, France to monitor the neutron 
scattering intensity in the spin-flip (SF) and non-spin- 
flip (NSF) channels at a number of high-symmetry points 
in the reciprocal lattice. Fig. 1 shows the SF intensity 
(normalized to the NSF intensity) at the (100) and (101) 
reciprocal-lattice vectors as a function of tempera- 
ture, with the neutron polarization P || Q. Temperature- 
independent contributions to the intensity are present 
at both positions. These contributions originate mostly 
from the leakage of NSF scattering from the correspond- 
ing nuclear Bragg reflections (which are not forbidden at 
these reciprocal- lattice vectors) , due to the usual limita- 
tions of the instrument. For Q = (101), an additional 
contribution to the SF intensity originating from elec- 
tronic magnetic scattering appears below Tn- We also 
observed magnetic intensity below Tn at (103), (105), 
(111), and (113), whereas none was observed at (100) 
[see Fig. 1] or (102), within experimental error. This 
indicates that the magnetic propagation vector is (001). 
The unit cell of 7-phase Na^OoO^ contains two planes 
of CoC>6 octahedra; a magnetic ordering vector of (001) 
therefore corresponds to AF ordering in the c-direction, 
combined with FM order within the ab- planes (Fig. 1, 
inset). The dependence of the SF intensity on the neu- 
tron spin direction at the sample confirms the inference 
from the uniform susceptibility of the c-axis orientation 
of the magnetic moment The (001) Bragg reflec- 

tion is unobservable because the spin orientation factor 
in the elastic neutron scattering cross-section vanishes for 
moments directed along c. By comparing the magnetic 
intensity at (101) with the nuclear intensity [2(| at (100), 
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FIG. 2: Scans at constant energy transfer hui taken along the 
(hhO) direction through (001) (for hui < 5 meV) or (003), 
at 1.6 K. The fits are to double Lorentzians. The red bar 
indicates an estimate of the FWHM of the instrumental reso- 
lution for the scan at hui — 5 meV. For clarity, the data have 
been offset vertically for the scans with hui < 15 meV (the 
corresponding vertical-axis zeroes are indicated at right) and 
the 15 meV data have been scaled by a factor of 2. Inset: 
constant-energy scans along (hhO) through (001) at energy 
transfer of 5 meV, at 1.6 K and 100 K. The fits are to double 
Gaussians. The axis units are the same as in the main panel. 



using the isotropic form factor, we extracted a value of 
0.13 ± 0.02 fi B per Co. 

The inelastic experiment was performed on the spec- 
trometer IN20 at the Institut Laue-Langevin, Greno- 
ble, France. We used an unpolarized configuration with 
a double- variable- focusing Si(lll) monochromator and 
a PG(002) analyzer with fixed vertical and horizontal 
focusing. In order to maximize the flux for the rel- 
atively small samples, which had volumes of 147 and 
162 mm 3 , no collimation was employed. Most of the 
data were taken either with a fixed final neutron wavevec- 
tor of kf=2.662 A -1 or a fixed incident wavevector of 
&i=4.1 A. -1 . The full width at half maximum (FWHM) 
of the energy resolution at zero energy transfer is approx- 
imately 1 meV in the former arrangement and 3 meV in 
the latter. The samples were each mounted in the (HHL) 
scattering plane and placed in a pumped 4 He cryostat. 

Low-temperature scans at constant energy transfer Tiuj 
through the ordering wavevector (001) along the in-plane 
(hh§) direction are shown in Fig. 2. Since the instrumen- 
tal resolution width is narrow relative to the peak widths, 
especially at large energy transfers, the data shown were 
fitted to unconvoluted Lorentzians. As a consequence of 
the large peak widths at higher excitation energies, in- 
terference from an optical phonon at hui = 20 meV, and 
the lower neutron flux at larger incident energies, it was 
not possible in this experiment to determine the disper- 
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FIG. 3: Constant-Q scans taken at different values of (0, 0, 3— 
0, with T = 1.6 K and k f = 2.662A" 1 . The curves are the 
results of convolutions of the full spin-wave dispersion (see 
Fig. 4) with the instrumental resolution ellipsoid; the blue 
bar indicates an estimate of the FWHM of the instrumental 
resolution for the scan with / = 0.25. 



sion relation up to the in-plane magnetic zone bound- 
ary. The inset shows similar scans taken with energy 
transfer fico = 5 meV at 1.6 K and 100 K, with coarser 
instrumental resolution (fcj = 4.lA _1 ). The integrated 
intensity decreases with increasing temperature, charac- 
teristic of a magnetic excitation. The separation between 
the two peaks is indistinguishable from zero at 100 K, 
which suggests that this signal arises from overdamped 
low-dimensional fluctuations. Comparison of this low- 
temperature scan through (001) with an equivalent one 
through (003) (not shown) demonstrates that the inten- 
sity decreases with increasing | Q | , as expected for a mag- 
netic excitation. No magnetic intensity was detected in 
similar scans through (002) and (004). 

For out-of-plane wavevectors, the magnon dispersion 
could be mapped out over the entire magnetic Brillouin 
zone. A set of constant-Q scans is shown in Fig. 3. 
The curves show results of a convolution of the full 3D 
magnon dispersion scattering with the instrumental res- 
olution, assuming a Lorentzian lineshape. Consideration 
of the fully convoluted dispersion was necessary for accu- 
rate determination of the peak positions because of the 
steep dispersion of the magnon branches in the orthog- 
onal (hhO) direction. In data taken at 50 K (~ 2.5Tn), 
the magnons are no longer present. 

The dispersion data in both directions are summarized 
in Fig. 4. Assuming nearest-neighbor interactions only, 
the spin Hamiltonian of an A-type antiferromagnet with 
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FIG. 4: a) Spin- wave dispersion along (hhO); the peak posi- 
tions are from Fig. 2. b) Dispersion along (00Z) , with peak 
positions from Fig. 3. All data were taken at 1.6 K. The dis- 
persion curves are the result of a fit to the model described 
in the text. 

localized spins can be written as 

H = J„ S ; -S j+ J ± £ Si-Sj-^S^, (1) 

(i,j) ihj) i 

in— plane .Lplane 

where Si is the spin-1/2 operator for the magnetic ion at 
lattice site i and the coupling constants Jm and Jj_ charac- 
terize the exchange interactions within the etfr-plane and 
between adjacent planes, respectively. The anisotropy 
constant D, the sign of which alternates from layer to 
layer, models the exchange anisotropy; it quantifies the 
tendency of the spins to align along the c-axis. The 
magnon dispersion can be calculated using the Holstcin- 
Primakoff formalism |20j| ; the resulting spin- wave disper- 
sion expression for q = (hkl) is 

= 2S yj {J ± (0) - [ J||(0) - Jiq)} + (D/2S)} 2 - {J ± (q)} 2 , 

(2) 

where j7j|(q) = 2J|| [cos(2tt/i) + cos(2ttA:) + cos(27r(ft, + k))] 
and i7j_(q) = 2J^cos(ttZ). We fitted the dispersion data 
in the (hhO)- and (OOZ)-directions simultaneously and 
extracted the following values: J|| = —4.5 ± 0.3 meV, 
J_L = 3.3 ± 0.3 meV, and \D\ = 0.05 ± 0.05 meV. Given 
the error in the peak position associated with the energy 
scan through the zone center at (003), and the large 
error in the fitted value of the anisotropy parameter, the 
question of whether there is an excitation gap at the AF 
zone center cannot be answered definitively by our data, 
and awaits further experiments using cold neutrons. Due 
to the dilution of spin-1/2 sites assumed for x — 0.82, the 
above calculation is only quantitatively accurate if the 
material is phase-separated into magnetically ordered 
regions with a dense network of spin-1/2 sites and 
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nonmagnetic regions. Another scenario, in which the 
spin-1/2 sites form an ordered superlattice, is discussed 
below. However, if the charge is uniformly distributed 
in the C0O2 layers or exhibits a small modulation in 
density, the calculation should be repeated using an 
itinerant model. 

Several aspects of our data are surprising. First, the 
spin-wave dispersion along the a6-plane is considerably 
steeper than that along the c-direction, but the magnon 
bandwidth is proportional to the number of nearest 
neighbors, which is six within the afe-plane, and only two 
along c. Hence, J|| and Jj_ are in fact comparable in 
magnitude. This relative isotropy is surprising in light 
of the two-dimensionality of the Na^OoC^ crystal struc- 
ture. In other layered magnets with comparable bond 
length anisotropics, such as YF^CuaOe, the magnitudes 
of the in-plane and out-of-plane exchange parameters dif- 
fer by orders of magnitude . It is also surprising that 
the Curie- Weiss temperature inferred from the fitted ex- 
change parameters in the context of a local-moment pic- 
ture, (6Jy + 2Jj_)/kB, is positive, whereas that extracted 
from the magnetic susceptibility measurements is nega- 
tive. 

The microscopic origin of these findings should be ad- 
dressed by theory. Since an unusually strong c-axis su- 
perexchange coupling through Na or an unusually weak 
nearest-neighbor in-plane exchange coupling would ex- 
plain the isotropy of the spin wave dispersions, a com- 
putation of these quantities is particularly desirable. A 
recent theoretical study has pointed out the relevance of 
longer-range exchange interactions along c • Another 
possible origin of our observation is an in-plane charge- 
ordered superstructure recently suggested pursuant to 
NMR 0, and optical conductivity experiments. 
The distribution of cobalt valence states has thus far not 
been elucidated directly. If all of the Co ions are in low- 
spin local-moment states (that is, S = 1/2 for Co 4+ and S 
— for Co 3+ ), then for x = 0.82 the spin lattice is dilute, 
with only 18% of the Co sites occupied by S = 1/2 spins. 
One possible arrangement of these spins, achievable ex- 
actly for x — 0.75, is that of a triangular lattice with 
lattice constant 2a. Since this is comparable to the out- 
of-plane Co-Co distance, isotropic spin wave dispersions 
would be a natural consequence. Such a spin ordering 
could be inferred from more complete measurements of 
the a6-plane magnon dispersion. A scenario in which the 
Co 3+ ions in the superstructure are in an intermediate- 
spin state with S = 1 and interact antiferromagnetically 
I Kill could help explain why the bulk DC susceptibility in 
Nao.82Co02 is AF, despite the fact that 6J11 + 2J_l > 0. 

If the cobalt valence is in fact distributed relatively 
uniformly, an itinerant-electron picture should be consid- 
ered instead. In this case, the magnetic ordering would 
correspond to a spin-density wave. 

Another interesting aspect of our results is the q- width 
of the (hhO) magnons, which increases with increasing |q| 



(Fig. 2). The large peak widths may reflect short-ranged 
magnetic correlations within the afe-planes. An alternate 
possibility is Landau damping by charged quasiparticles. 
In a charge-ordered scenario with antiferromagnetically 
interacting, intermediate-spin Co 3+ ions arranged on a 
geometrically frustrated Kagome lattice 0], the broad- 
ening may also arise from an admixture of excitations 
from the disordered array of Co 3+ spins. 

In conclusion, our measurements demonstrate that 
Nao.82Co02 exhibits spin fluctuations characteristic of 
low-tempcrature 3D AF ordering of the A type, with 
magnetic exchange constants much less anisotropic than 
expected, given the layered crystal structure. The 
antiferromagnetism coexists with metallic conductivity. 
Unanswered questions remain regarding the degree to 
which the spins have localized or itinerant character, pos- 
sible superstructure in the former case, and the micro- 
scopic character of the exchange couplings. Answers to 
these questions may also shed light on the origin of the 
unusual thermal properties of Naa,Co02 and on the su- 
perconducting state in the hydrated analogue. 
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